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We analyze the 1 + 1 dimensional Nambu-Jona-Lasinio model non-perturbatively. We 
study non-trivial saddle points of the effective action in which the composite fields 
oix) = (V''!/') and ■k(x) = {'ipi^y^ip) form static space dependent configurations. These con- 
figurations may be viewed as one dimensional chiral bags that trap the original fermions 
("quarks") into stable extended entities ("hadrons"). We provide explicit expressions 
for the profiles of some of these objects and calculate their masses. Our analysis of 
these saddle points, and in particular, the proof that the cr{x), tt{x) condensations must 
give rise to a reflectionless Dirac operator, appear to us simpler and more direct than 
the calculations previously done by Shei, using the inverse scattering method following 
Dashen, Hasslacher, and Neveu. 

1. Introduction 

In this talk we describe a novel methodl for studying the non-oerturbative spectrum 
of the 1-1-1 dimensional Nambu-Jona-Lasinio (NJL) model0Q 
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The action (1.1) describes N self interacting massless Dirac fermions V'a (a = 
1, . . . , N), and we study it in the limit N ^ oo holding Ng'^ finite (the large N 
limit). This model is interesting because it shares with QCD some of its important 
low energy attributcsa, and at the same time its non-trivial dynamics is under con- 



trol. Indeed, (1.1) is invariant under SU{N)f (g) U{1) (g) U{1)a- It is asymptotically 
free and exhibits dynamical mass generation (for iV > 2) due to its infra-red insta- 
bilities. These instabilities may also polarize the vacuum dnhomogenously, giving 
rise to chiral solitonaJ'El that are reminiscent of chiral baggj'Ll which model hadron 
formation. These chiral solitons and the small fermionic fluctuations around them, 



are the main subject of this talk. We rewrite (1.1) as 



j (fx [i$ - (a + z^75)] V' - ^ (fr" + ^')| 



(1.2) 



*Talk delivered by J.Feinberg at the Workshop on Low Dimensional Field Theory at Telluride, CO 
(August 1996). To appear in the conference proceedings. 
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where cr(x), ■k{x) are the scalar and pseudoscalar auxiliary fields, respectively, Which 
are both of mass dimension 1. These fields are singlets under SU{N)f ® U{1), 
but transform as a vector under the axial transformation, namely a + 47577 — > 
g-2j75/3j^^ _l_ j^gjr) ^ Gaussian integration over the fermions in ( |1.2| ) leads to the 
partition function Z = J DaDn exp{iS'e//[cr, tt]} where the bare effective action is 



1 



(fx (cr^ + TT^) - iN Tr In [i0 - {a + ^75)] 



(1.3) 



and the trace is taken over both functional and Dirac indices. 



The non-perturbative vacuum of (1.3) is determined, in the large N limitH, by the 



simplest large N saddle points of (1.3) where a and tt develop spacetime independent 
expectation values. These saddle points are extrema of the effective potential Ve//, 
which as a result of chiral symmetry depends only on the combination = cr^ +7r^. 
It is minimized at p = mA e '''a'-'w ^ where m is tl^ dynamical mass of fermions. 



This mass is determined by the (bare) gap equatiorn 



m + iNg tr 



(1.4) 



were A is an ultraviolet cutoff. The mass m must be a renormalization group 



invariant. Thus, the model is asymptotically free. The vacuum manifold of (1.2) is 
therefore a circle p = m in the u, tt plane, and the equivalent vacua are parametrized 
by the chiral angle = arctan^. Therefore, small fluctuations of the Dirac fields 
around the vacuum manifold develop dynamical-'-chiral mass to exp(i6'75). § 
Static, space dependent solutions of the saddle point equations of (pT^) 



SS. 



Sa {x, t) 
Sir {x, t) 



a {x, t) 

TT [x, t) 



iN tr 



N tr 



{x, ti- 
ls {x,t\ 



1 



id — {(J + 47775) 
1 



W — (cr + 17775) 



|x,t) 

\x,t) 







(1.5) 



are the chiral solitons (at rest) mentioned above. The NJL model, with its contin- 
uous symmetry, does not have topologically stable soliton solutions. The solitons 
arising in the NJL model can only be stabilized by binding fermions and releasing 
binding energyliEI. This observation is clarified further by comparing the NJL model 



to the Gross- NeveuH (GN) model, namely, (1.1) with the 75 term deleted. The GN 
model possesses only a discrete symmetry, ijj —f "f5ip, rather than the continuous 
symmetry of the NJL model. This discrete symmetry is dynamically broken by 
the doubly degenerate non-perturbative vacuum, and thus there is a topologicaly 

tProm this point on flavor indices are suppressed. 

•I-Note that the axial U{1} symmetry protects the fermions from developing a mass term to any 
order in perturbation theory. 

§Note in pmssing that the massless fluctuations of 6 along the vacuum manifold decouple from the 
spectriHp □ so that the axial U{1) symmetry does not break dynamically in this two dimensional 
model 
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stable kink solutiontlUlSO, the so-called Callan-Coleman-Gross-Zee (CCGZ) kink, 
that interpolates between the two vacua. The kink binds any number n < TV of 
fermions in its single zero energy bound state, without affecting its mass. Its sta- 
bility is guaranteed by topology already. In contrast, the stability of the extended 
objects arising in the NJL model is not due to topology, but to dynamics. 

Shei B hasjound certain static solutions of (1.5) by applying the inverse scat- 
tering methodliJ, following a similaranalysis by Dashen et al£3. Recently, one of 
us developed an alternative methodll3'll3,based on the Gel'fand-Dikii identityll3, to 
investigate spectra of such low dimensional field theories. We feel that this method 
has certain advantages over the inverse scattering method for analysing this type of 
problems, as it bypasses some of its technical steps. This talk is based on a paper i 
where we have applied this new method to study (LI). In particular, our proof that 
the Dirac operator in the background of the static solutions of (^]^) is reflectionless 
is extremely simple, and does not require invoking inverse scattering techniques. It 
is actually valid for finite N. It is worth mentioning at thisjioint that the NJL 
model (1.1) is completely integrable for any number of flavorsE-3 and thus its QS^-ct 
spectrum and S matrix are known in detail_The large N spectrum we describeBtli 
consistent with the relevant exact resuts in t3. However, the powerful methods of ll, 
are inherently limited to 1 + 1 dimensional integrable models, whereas the analysisi 
we describe here is potentially applicable to 1 -I- 1 dimensional non-integrable mod- 
els as well aa to analysing inhomogeneous symmetric field configurations in higher 
dimensionsta. 
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2. Absence of Reflections in the Dirac Operator With Static Back- 
ground Fields 

As was explained in the introduction, we need a manageable form of the diagonal 
resolvent of the Dirac operator D = {(j{x) + iTi{x)^^) in a given background of 
static field configurations (j{x) and t:{x). The extremum condition on S^g relates 
this resolvent, which in principle is a complicated and generally unknown func- 
tional of <j{x), 7r(a;) and of their derivatives, to ct{x) and 7r(a;) themselves. This 
complicated relation is the source of all difficulties that arise in any attempt to 
solve the model under consideration. It turns out, however, that basic field theo- 
retic considerations, that are unrelated to the extremum condition, imply that D 
must be refiectionless. This spectral property of D sets rather powerful restrictions 
on the static background fields cr(a;) and ■n{x) which are allowed dynamically. In 
the next section we show how this special property of D allows us to write explicit 
expressions for the resolvent in some restrictive cases, that are interesting enough 
from a physical point of view. Inverting D has nothing to do with the large N 
approximation, and consequently our results in this section are valid for any value 
of TV. 

The overall energy contained in any relevant static cr, tt configuration must be 
finite. Therefore these fields must approach constant vacuum asymptotic values. 



namely, points on the circle 



+ ±oo 



with vanishing derivatives. We use the 



axial U{1) symmetry to fix the coordinates in the cr, tt plane such that a{—oo) = m 
and 7r(— c») = 0. Then a{oo) — mcosO and 7r(oo) — msmO, where 9 is the chiral 
alignment angle of the vacuum at a; = +oo relative to the vacuum at x = — oo. We 
use the Majorana representation = cr2 , 7^ = ia^ and 7^ — —7^7^ = cri for 7 
matrices. In this representation D becomes 



D = 



— dx — cr — itiJ — ITT 



iuj — ITT dx — u 



(2.1) 



We invert (2T) by solving 



—dx — <y{x) —iuj — iTT{x) 
iu) — i-K{x) dx — o'ix) 



a{x,y) b{x,y) 
c{x,y) d{x,y) 



-il6ix-y) (2.2) 



for the Green's function of (2.1) in a given background a{x), n(x). By dimensional 



analysis, we see that the quantities a, 6, c and d are dimensionles. We analyzed (2.2) 
in 13. The quantities a[x,y),c{x,y) and d(x,y) in (2.2) may all be inferred from 
b{x, y), through symmetries of ( |2.lD . The quantity b{x, y) is the Green's function of 
the one dimensional Sturm-Liouville operator 



-dx 



dxb{x,y) 

LU + Tt(x) 



(j{xf + 'K{xf - cr'(x) -ijj"^ + 



a{x)T:' {x) 

LO + TT^x) 



b{x,y) 
o + 7r(a;) 



= ^x-y). 
(2.3) 
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Analysis of the asymptotic behavior of ( |2.2| ) and (2.3) turned out to be quite 
valuable. The physical asymptotic boundary conditions imposed on cr, tt are such 
that the potential term in (2.3) asymptotically a constant, — uj'^. We therefore 



deduced in i that the asymptotic behavior of the diagonal resolvent of (2.1) is 

1 + i? (fc) e^^*^ 1^1 



iD- 



X — »±oo 



2fc 



[«757r(a;) - a{x) 



R{k)e 



2ik \x\ 



■ 7 sgn X 



(2.4) 



where k = \/ uP- — vn? and i?(fc) is the reflection coefficient of the Sturm-Liouville 
operator in (2.3). Consider now the expectation value of fermionic vector current^ 
operator in the static cr(a;), 7r(a;) background 



(cr(x), 7r(a;)|j^|f7(a;), 7r(x)) 



27r 



tr 7*^ 



iD- 



(2.5) 



We thus find from ( ]2.4D that the asymptotic behavior of the current matrix elements 
is such that {(j{x), 7r(x)|j''|tT(x), 7r(a;)) — > 0, but 



(cr(x), 7r(a;)|/|cr(x), 7r(a;)) 



/ gi?(fc)e2^'=l^l sgn: 



(2.6) 



where we used the fact that J §^^f{k) — because k (w) is an even function of uj. 
Thus, an arbitrary static background cr(x) , 7r(a;) induces space dependent fermion 
currents that do not decay fast enough as x — > ±oo, unless R{k) = 0. Clearly, we 
cannot have such currents in our static problem and we conclude that as far as the 



field theory (1^) is concerned, the fields a{x) ,Tr{x) must be such that the Sturm- 
Liouville operators in (2^) and therefore the Dirac operator (2T) are reflectionless. 

One may also write equations for (cr(x), TT{x)\^p'ip\a{x), tt{x)) and 
{cr{x), 7r(a;)|iV'75V'|o'(a:), n{x)) that are analogous to (2^). For reflectionless tr, tt 
backgrounds, these equations boil down at the saddle point (where the expectation 
value of the scalar density is — cr(a;)/.g^, and that of the pseudoscalar density is equal 
to —Ti{x)/g'^) into 

1 ^, , 1 

(2.7) 



Ng' 



1 



duj 



27r \/llP- — rn? 



which is simply the gap equation (1.4), and thus holds to begin with 
The absence of reflections in (2.1 



L) emerges here from basic principles of field 
theory, namely, that an initially static (7(a;), ^{x) will not become time dependnt 
due to backreaction. Absence of reflections is thusL valid beyond the large iV saddle 
point condition from which it was deduced in I13i3. For reflectionless backgrounds 
( ^■4| ) simplifies to 

1 



iD- 



■>±oo 2V^ 



\i^^T^{x) - a{x) - 7°^] 



(2., 



■^In the following it is enough to discuss only the vector current, because the axial current jj" 
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This expression has cuts in the complex uj plane corresponding to scattering states 
of fermions of mass to. These cuts must obviously persist away from the asymptotic 
region, and we make use of this fact in the next section. 



3. The Diagonal Resolvent for a Fixed Number of Bound States 



The requirement that the static Dirac operator (2T) be reflectionless is by itself 
quite restrictive, but in order to actually derive explicit expressions for the resolvent 
in terms of cr{x) , 7r(a;) and their derivatives we supplement it by assuming in addition 
that a{x) and t:(x) are such that the spectrum of contains a prescribed number 
of bound states. In the following we concentrate on the diagonal resolvent B{x) = 

x) of (|2^ 
Dikii identit\ 




Being the diagonal resolvent of (|2.3| ), B{x) satisfiesH the Gel'fand 



1 



Lo + 7r(a;) 



+ 7r(a;) 
2B{x) 



d:,B{x) 
u! + Tr{x) 



B{x) 



UJ + tt{x) 



a{xf + n{xf - a'{x) - u;^ + ^MZlM 

• ■ ' t:(x) 



to ■ 



Tl{x)Y 



■ dx 



a{xf + t:(xY - a'{x) + 



cr(x)7r'(a:) 



it{x) 



= 0. 



(3.1) 



If we were able to solve (3.1) for B{x) in a closed form for any static configuration of 
(j{x), 7r(a;), we would then be able to express {x \iD~^ \x ) in terms of the latter fields 
and their derivatives, and therefore to integrate (|l.5|) back to find an expression for 



the effective action (1.3) explicitly in terms of (j{x) and t:{x). Invoking at that point 
Lorentz invariance of (1.3) we would then actually be able to write down the full 
effective action for space-time dependent a and tt. Note moreover that in principle 
such a procedure would yield an exact expression for the effective action, regardless 
of what N is. Unfortunately, deriving such an expression for B{x) in general is a 
difficult task. It becomes manageable only for the special cr, tt backgrounds specified 
above. 

We now sketch our treatments of a reflectionless background that supports a 
single bound state. In such a case, B{x^ lo) has to have a simple pole on the real uj 
axis, say at w = wi, where clearly, \loi\ < m. There is also the continuum cut as in 
(2.8). The product of these two factors has dimension —2 in mass units. Any other 



singularity B(x,uj) may have in the complex lu plane must involve x dependence as 
well, through the combination exp(z-\/w^ — m? x). However, the requirement that 
the Dirac operator be reflectionless, rules this possibility out. The pole and cut then 
exhausts all allowed singularities of B(x,lu) in the complex uj plane. In addition, 
B{x,uj) is dimensionless, and therefore has to be of the form 

B2ix)uj^ + Bi{x)uj + Baix) 



B{x,Lu) 



V TO^ — LlI^{llJ — LlJi) 



(3.2) 



where the dimension of the unknown function B^ix), (fc = 0, 1, 2) is 2 — k. These 
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functions are combinations of a{x),TT(x) and their derivatives. Substituting (3.2) 



into (3.1) results -^in a polynomial of degree six in uj, with x dependnet coefficients 
that has to vanish identically. In this way they form an over-determined set of seven 
differential equations in the five functions B2, Bq, a{x) and tt(x). However, this 
over-determined system has a solutionEI satisfying the boundary conditions (2.8), 



„, ^ UJ + Tr(x) a'^(x) + Tr'^(x) ~ a (x) ~ m'^ 
B{x) = ^ ' + — ^ ' ^ ' where (3.3) 



2V^ 



4 (tj — ciJi) \J 



(j{x) = TO — 



msin^tan' 



1 -I- exp [2tJitan| • (x — a;o)] 



and 



(3.4) 



7r(x) 



\(y{x) — m] cot— . 



(3.5) 



This linear relation between a and tt is not surprising, because tr, tt are the two 
components of an axial vector. Note that the boundary conditions at a; — s- -l-oo 
require 



Wi tan- < . 



(3.6) 



Our results (3.4) for o{x) and 7r(x) agree with those of B. They have the profile of 
an extended object, a lump or a chiral "bag" , of size of the order cot|/a;i centered 



around an arbitrary point xq. Note that the profiles in (3.4) satisfy 



p^{x) = CT^(x) + Ti^ix) 



m 



TO^ sin^ (6'/2)sech2 



Witan- • {x — Xq) 



(3.7) 



Thus, as expected by construction, this configuration interpolates between two dif- 
ferent vacua at a; = ±cxd. As x increases from —00, the vacuum configuration be- 
comes distorted. The distortion reaches its maximum at the location of the "bag" , 



where 



p^{xo) — TO^ sin {0/2) and then relaxes back into the other vacuum 



state at a; = cxD. The arbitrariness of xq is, of course, a manifestation of translational 
invariance. 

In i we also have partial results for the case of two bound states, at energies 
LUi, UJ2- In particular we found that y{x) — (T{x)^m satisfies the differential equation 



2X{LJi+ U2) [4m + 2(1 + A^)^^ - y"] 

+84 4 (to^ + UJ1UJ2) y + 6my^ + 2{X^ + l)y^- y"} = 



(3.8) 



where A — — cot(0/2), and 7r(x) is given by (3.5). Note that if we set uji + 0J2 = 
and Tr{x) — the spectrum becomes invariant under lo — > —u, and we obtain the 
equation appropriate to the Gross-Neveu model. 



^Nol£that because of the linearity and homogeneity of (3.1), the purely lu dependent denominator 
of (pj) with its explicit dependence on the bound state energies factors out. 
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4. The Saddle Point Conditions 



Derivation of the explicit expressions of <j{x) and 7r(x) does not involve tlie saddle 
point equations (|1.5|). Rather, it tells us independently of the large N approximation 



that (t{x) and ^{x) must have the form given in (3.4) in order for the associated 
Dirac operator to be reflectionless and to have a single bound state at a prescribed 



energy loi in addition to scattering states. Thus, for the solution (3.4) we have yet 



to determine the values of uji and 9 allowed by the saddle point condition (1.5). It is 



this dynamical feature that we can analyse only in the large N limit. Note that the 
from (3.3) we can reconstruct the other three entries of {x \ —iD^^\x) . Substituting 



this resolvent, with cr, tt given by (3^) into (|l.5| ) we obtain a quantization condition 
for wi , in the form of coupled dispersion integrals that have to vanish. A nice feature 
of these integrals is that their potentialy ultra-violet divergent parts are equal to 
the difference of the two sides of the gap equation (^), and are therefore free of 
such divergences. Some lo poles in these integrals are at x dependent locations. The 
quantization condition on wi cannot be x dependent, so the residues of these poles 



must vanish. This happens provided the parameters in (3.4) satisfy 



^ — m^cos^(-) 



(4.1) 



This relation actually leaves the only free parameter in the problem with respect 



to which we have to extremize the action. The condition (3.6) then picks out one 
branch of (4J). Assuming that the Dirac sea is completely filled, and that the 
bound state traps n < N fermions, the saddle point condition becomes equivalent!! 
to the requirement that Seff evaluated at the appropriate cr, tt configurations, be 
stationary as a function of 9, namely. 



1 dS _m 



Only the critical value at 



n 
N 

27m 



I sm- . 
27r / 2 



(4.2) 



(4.3) 



corresponds to extended objects. The other zeros of (4.2) at 9k — 27rfc, A: g Z do 



not correspond to such objects at all as should be clear from (3.4). We therefore 



discard them, and concentrate on the extremum given by (4.3). The bound state 
energy is therefore lui = mcos (7^ ). Equation ( |4.3| ) asserts that the relative chiral 
rotation of the vacua at ±00 is proportional to the number of fermions trapped in 
the "bag". Our result ( |4.3| ) is consistent with the fact that the fermion number 
current in a soliton background can be determined in some cases by topological 
considerationsEa. Note from (4.3), that in the large N limit, 9 (and therefore uji) 
take on non-trivial values only when the number of the trapped fermions scales as 
a finite fracion of N. 

As we already mentioned in the introduction, "bags" formed in the NJL model 
are not stable because of topology. They are stabilized by releasing binding energy 
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of the fermions trapped in them. To see this more exphcitly, we calculate now the 



mass of the "bag" corresponding to (B.4) and (4.3). Integrating (4.2) with respect 
to 6 we find 



-1 



NTm 



3(0) 



n 
N 



, cos- 
27r / 2 



1 . 
— sm— 

TT 2 



(4.4) 



Note that ( [4.4[) is not manifestly periodic in 9 because the Pauli exclusion principle 
limits 9 to be between and 2tt. The mass of a "bag" containing n fermions in 
a single bound state is given by the energy E{9) = —S{9)/T evaluated at the 
appropriate chiral angle (O). We thus find that this mass is simply 



Nm TTTl 

sm— — 

TT N 



(4.5) 



in accordance with EHIII. It is easy to check that (^ 
< < 1. These "bags" are stable because 



T) is a minimum of E{9) for 



. TT (ni + 712) . TTUl 

sm < sm 

N N 



sm- 



7rn2 



(4.6) 



for rti + n2 less than N, such that a "bag" with ni + 712 fermions cannot decay into 
two "bags" each containing a lower number of fermions. 

Entrapment of a small number of fermions cannot distort the homogeneous 
vacuum considerably, so we expect that M„ will be roughly the mass of n free 
massive fermions for n << iV as (15) indeed shows. As a matter of fact we used 



this expectation to fix the integration constant in (4.4). However, as the number 
of fermions trapped in the "bag" approaches N, Mn vanishes and the fermions 
release practically all their rest mass Nm as binding energy, to achieve maximum 
stabilityllj. Note from (|4.3| ), that the soliton twists all the way around as the number 
of fermions approaches N. 

We conclude this talk making some remarks on the case of two bound states. Ini 
we found that the bound state energies must be located atuoi — m cos (^^) , t^2 ~ 
rn cos (^^) which are identical in form to single bound state energy levels. From 
the general considerations of Ej we expect that the chiral angle 9 will be propor- 



tional to the total number of fermions trapped by the "bag", so (4.3) must read 



27r(ni+n2) 
N 



Clearly, the mass of the "bag" will depend on 711 and n2 



seperately, and not only on their sum through 9. 
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